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ABSTRACT 

Physical picture of the emission mechanisms operating in the X-ray binaries was put under question by the 
simultaneous optical/X-ray observations with high time resolution. The light curves of the two energy bands 
appeared to be connected and the cross-correlation functions observed in three black hole binaries exhibited 
a complicated shape. They show a dip of the optical emission a few seconds before the X-ray peak and the 
optical flare just after the X-ray peak. This behavior could not be explained in terms of standard optical emission 
candidates (e.g., emission from the cold accretion disk or a jet). We propose a novel model, which explains 
the broadband optical to the X-ray spectra and the variability properties. We suggest that the optical emission 
consists of two components: synchrotron radiation from the non-thermal electrons in the hot accretion flow and 
the emission produced by reprocessing of the X-rays in the outer part of the accretion disk. The first component 
is anti-correlated with the X-rays, while the second one is correlated, but delayed and smeared relative to the 
X-rays. The interplay of the components explains the complex shape of the cross-correlation function, the 
features in the optical power spectral density as well as the time lags. 

Subject headings: accretion, accretion disks — black hole physics — methods: numerical — X-rays: binaries 



1. INTRODUCTION 

The physical processes giving rise to the broadband radio to 
X-ray spectra of accreting black holes are still under debate. 
There is a general consensus that the radio emission is associ- 
ated with the jet and the X-rays are produced in the vicinity of 
the compact object. However, the origin of the IR/optical/UV 
emission from the black holes in low-mass X-ray binaries 
(LMXBs) is less certain. The contribution from the compan- 
ion star in such systems is usually too faint, and the optical 
spectra are most likely connected to the accretion process onto 
the compact object and the X-ray radiation. To understand the 
origin of the optical emission, timing analyses of simultane- 
ous optical and X-ray observations were performed. S uch ob- 
servati ons were for the first time carried out by Motch et al.l 
( [1981 for the black-hole binary (BHB) GX 339^. Although 
the duration of the observations was too short for any confi- 
dent conclusion, the computed optical/X-ray cross-correlation 
function (CCF) revealed a complicated structure with a dip in 
the optical light curve, preceding the X-ray peak (the so-called 
precognition dip), together with an optical peak lagging the 
X-rays. Recently, CCFs were obtained from the much longer 
duration simultaneous observations in three LMXBs: XTE 
J1118H-480 (Kanbachetal. 2001; Hynes et al. 2003), Swif t 
J1753.5-0127 (Durant et al. 2008, 20111: iHvnes etal]l2009h 
and GX 339^ (Gandhi et al. 2008^), all manifesting similar 
structure. The observed behavior cannot be explained by a 
simple model with the optical radiation being pro duced by 
the reprocessed X-ray emission (iHvnes et al.ll2003l) . This hy- 
pothesis also contradicts the fact that the observed autocorre- 
lation function ( ACF) of the optical radiat i on is narrower than 
that of the X-rays (Kanbach et al. 200 Ij iHynes et all 120031: 
iDurant etal., .2009: G andhi et al. 2010) . On the other hand, 
study of the dependence of the CCF on the timescale of fluc- 
tuations in the light curve r evealed that its overall shape re- 
mains similar, but rescaled (iMalzac et al.ll2003l: iGandhi et alj 
|20T3), suggesting that the emission in the two energy bands 



is intrinsically connected. 

A number of possible mechanisms, producing such struc- 
ture of the CCF , were proposed (see the discussion in 
iDurant et al.ll201 1 and references therein). So far, the d etailed 
calcul ations were made only in the model of Malzac et alJ 
(l2004l) . where both X -ray and optical emission are powered 
by the same magnetic energy reservoir The model qualita- 
tively describes the CCF and ACF of XTE Jl 1 18H-480, how- 
ever, with the new data on Swift J1753. 5-0127, the applica- 
bility of the model was put under question. It is also doubtful, 
that such an energetic magnetized zone is consistent with the 
accretion theories. 

In this Letter, we suggest that the optical emission is a com- 
position of two components. The first one is coming from 
the synchrotron-emitting particles in the inner hot flow and is 
anti-correlated with the X-ray emission, which is produced 
by Comptonization of the synchrotron radiation. The sec- 
ond one is originating from reprocessin g of the X-rays by 
the outer pa rts of the col d accretion disk (iDubus et al.llI999t 
iGierlinskiet al. 2009; Chi ang et al.ll2010l) . This component is 
(positively) correlated and delayed relative to the X-rays. The 
presence of two components explains both the precognition 
dip and the delayed peak of the CCF. The model is capable of 
explaining the entire optical to the X-ray spectrum of BHBs 
as well as their timing properties. 

2. SPECTRAL MODEL 

The X-ray spectra of BHBs in their hard state are 
well described by C omptonizati on on therrnal e lectrons 
(e.g., iPoutanenI II998I; [Zdziarski & GierUnskil 120041) . how- 
ever power-law tails detected in a number of sources 
(^dziai-ski et al. 2001; McConnell et al. 2002) suggest the 
presence of non-thermal particle population in addition t o 
the mostly thermal distribution (iPoutanen & CoppH Il998h . 
These electrons contribute both to the high-energy part of 
the spectrum via Compton scattering and to the low-energy 
part by synchrotron radiation. Even a tiny fraction of non- 
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Figure 1. Simulated (a) photon spectra and (b) electron distiibutions for 
the SSC model for hybrid plasmas. Solid lines correspond to the fiducial 
parameter set (see the text), an increase of the accretion rate by 20% hardens 
the spectrum (blue dashed line) and a decrease of the accretion rate by 20% 
softens the specfium (red dotted line). These spectra are associated with the 
inner hot flow and an additional contribution of the irradiated accretion disk 
is not shown. 

thermal electrons dramatically increases the synchrotron lu- 
minosity (Wardzinski & Zdziarski 2001) and thus makes it a 
plausible source of seed photons for Comptonization. The 
X-ray spectra of the hard state BHBs (and supermassive 
black holes) can be modeled in terms of synchrotron self- 
Compton (SSC) mechanism in hybrid (thermal plu s non- 
therm al) pl asmas (Vurm & Poutanen 2008; Poutanen & VurmI 
l2009t iMaizac & Belmont ,20091: iVeledina et al.i i201 ID . The 
high-energy part of the spectrum is dominated by the Comp- 
tonized photons, and the low-energy part is determined by 
the synchrotron emission, which can extend down to the 
IR/optical energy bands depending on the size of the emis- 
sion region. 

For simplicity, we consider a spherical emission region ho- 
mogeneously filled with photons and electrons. The region 
corresponds to a hot flow in the vicinity of the blac k hole 
(see reviews in iNaravan etallfTOQl lYuanet al.ll2007h . wifli 
the size limited by the inner radius of the (truncated) accre- 
tion disk. We consider the energy input in the form of injected 
electrons with a power-law spectrum dN^/ {dt dj) oc 7^^""'^ 
extending between the Lorentz factors 7inin and 7max- The 
injection might result from magnetic reconnection or shock 
acceleration. The main mechanisms responsible for particle 
cooling and thermalization and formation of the spectra are 
cyclo-synchrotron emission and absorption, Compton scatter- 
ing and electron-electron Coulomb collisions. To find self- 



consistent photon spectra and electron distributions, we solve 
a set of kinetic equations for photons, elect rons, and positrons 
using the numerical code developed by IVurm & PoutanerJ 
(|2009l). 

A good description of the broadband spectra of hard-state 
BHBs can be achieved with the following parameters (see 
iPoutanen & Vurmll2009l) : the Thomson optical depth r=1.0, 
the magnetic field B = 3 x 10^ G and the total luminosity of 
the region L = 10'^'' erg s~^. The parameters of the injection 
function are Finj = 3, 7inin — 1-0, and 7inax = 10"^. The size 
of the emission region R is assumed to be 30 Schwarzschild 
radii, which corresponds to 9 x 10^ cm for a IOMq black hole. 
For this set of parameters, the synchrotron self-absorption fre- 
quency falls in the optical energy band (1 — 10 eV). The re- 
sulting spectrum and the electron distribution are shown in 
Figure [T](solid line). 

The broadband spectra are expected to vary due to the 
changes of the mass accretion rate m. We assume that the 
luminosity varies as L oc m, the optical depth as t (x m, and 
the size of the hot flow (inne r cold disk radius) as i? oc rh^"'/'^ 
dRozaiiska & Czemvl 120001) . We assume that the magnetic 
field and the injection index Finj are constant. With an in- 
crease of the accretion rate, the X-ray luminosity increases, 
whereas the optical radiation drops due to the increased syn- 
chrotron self-absorption and vice versa (dashed and dotted 
lines in Figure[T]i. The resulting X-ray spectral behavior with 
hardening during the X-ray flares and softeni ng during the 
X-ray dips is consistent with the observations (iGandhi et alJ 
I2OIO). The decrease of the synchrotron radiation with an in- 
crease of the mass accretion rate was also foun d to take place 
in the advective disk models dYuan et alj|2005h . Thus, the op- 
tical synchrotron luminosity is expected to anti-correlate with 
the X-rays. 

3. TIMING MODEL 

3.1. Formalism 

The X-ray and the optical light curves can be represented 
as a sum of a constant (mean) and variable components: 
x{t) = x[l + 5x{t)] and o{t) = d[l + So{t)], where Sx and 
So are the relative deviations from the mean. Two compo- 
nents contribute to the variations of the optical radiation: the 
synchrotron emission from the hot flow Ss and the radiation 
coming from reprocessing of the X-rays in the accretion disk 
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The anti-correlation between synchrotron and the X-rays (see 
Section|2]i can be represented as 

Ss{t) = -l^fa(t), (2) 
rmSx 

where rms denotes the fractional root mean square amplitude 
of variability of the corresponding components^] The repro- 
cessed radiation from the disk is a convolution of the X-ray 
light curve with the disk response function r{t): 



Sd{t) 



r{t - t')5x{t')dt' 



(3) 



Assuming the entire disk luminosity in the optical band is due 
to irradiation and the response is a (5-function, the fractional 

' We assume here that the decrease of the synchi'otron luminosity is simul- 
taneous with the increase of the X-rays, neglecting possible delays (see the 
discussion in Section|4). 
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Figure 2. Results of simulations for the model with the exponential disk response function with t\ =0.1 and T2 = 1.0, and r^s = 1.3. (a) The PSDs of 
the X-ray (double-broken PSD, thick solid line) and optical hght curves (thin solid Hne). Thi'ee terms contributing to the optical PSD (see Equation (sj) are also 
shown: synchrotron (dashed line), disk (dotted line), and the cross term (positive contribution is shown by the dot-dashed line and the negative contribution by 
the long-dashed line), (b) The CCFs, (c) the phase lags, and (d) the time lags coiTesponding to the synchi'otron/X-rays (rds = 0, dashed line), reprocessed disk 
component/X-rays (r^g — > oo, dotted line), and the combined optical/X-rays (solid line). 



rms of the disk would be equal to that of the X-rays. Total 
optical light curve can be represented as 



boit) 



rmSx o 



Sx{t) + Tds / 6x{t')r{t ~ t')dt' 



(4) 

where we introduced the ratio of the absolute contributions of 
the disk and the synchrotron to the variable optical component 
(with the assumption of a J-function response): 



rds = --■ (5) 

rmSs s 

As follows from Equation the optical/X-ray CCF contains 
two terms: one arising from synchrotron (it is just an ACF of 
the X-rays with a negative sign) and another from the disk. 

The Fourier transformsQ for the optical light curves are the 
sum of the synchrotron S{f) oc —X{f) and the disk D{f) — 
X{f)R{f) transforms, and according to Equation (|4) can be 
represented as 

0{f) = S{f) + D{f) « X{f) [-1 + r^sRif)] ■ (6) 

^ The Fourier transforms corresponding to Sx(t), So(t), Ss(t), 8d{t), and 
r{t) are denoted as X{f), 0{f), S{f), D{f), and R{f), respectively, with 
/ being the Fourier frequency. 



The optical/X-ray cross-spectrum is proportional to the X-ray 
power spectral density (PSD, Fx) 

C{f) = X*{f)0{f) cx Pxif) [-1 + rdsi?(/)] . (7) 

Its phase is the phase lag A0(/) and the time lag is At(/) = 
A0(/)/27r/. The optical PSD consists of three terms: the 
synchrotron term, the disk term, and a cross term 

Poif) « Pxif) {l + rl \Rif f - 2rdsRe [i?(/)]} . (8) 

The model is fully determined by the rms in the X-ray and 
optical bands (rmSx and rmSo), the ratio r^s, the disk response 
function r{t) and the shape of the X-ray PSDI3 

3.2. Examples 

For an illustration, we choose the double-broken power-law 
X-ray PSD, Pxif) oc /°, f~^ and f'^ with the breaks at 0.1 
and 1.0 Hz. This shape is typical for BHBs in their hard state 
(Gilfanov et al.' 1999; Rev nivtsev et all 120011: lAxelsson etal] 
2005; Gandhi et al. 2008J. We take rmSx = 1.0 and rmSo = 
0.3 and consider a simple exponential response function 



exp[~{t ~ Ti)/t2] /t2, t^Tl, 
0, t<Tl, 



(9) 



We choose such a PSD normalization so that the integral over pos- 
itive fre quencies gives the square of relative rms of the hght curve 
jMivamot o & Kita motclfT^ / Px{f)df = rms^. 
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Figures, (a) X-ray and optical PSDs, (b) optical/X-ray CCFs, (c) optical/X-ray phase lags, and (d) time lags for various r^g. The results for r^g =0.7, 1.0, 
1.3, and 2.5 ai'e shown by dashed, dotted, solid and long-dashed lines, respectively. The X-ray PSD is shown with the thick solid Hne. Other parameters ai'e as in 
Figure 12] 



where t\ corresponds to the delay time of the disk response 
and T2 is the response width. The corresponding Fourier im- 
age is exp(ixi)/(l — \x-2), where x-^ — 2TTfTk. We 
take Ti = 0.1 s and T2 = 1 s. The results of simulations are 
shown in Figure |2] 

Contributions of different terms and the total PSDs are 
shown in Figure |2a). The synchrotron term (the first term 
in Equation dHJ) has the same shape as the X-ray PSD, the 
disk (the second term in Equation (|8]l) acts as a low-pass fil- 
ter, and the cross term (the last term in Equation (|8}) changes 
sign, giving positive contribution to the high-frequency part 
and negative contribution to the low-frequency part. The com- 
bined contribution of the two (synchrotron and disk) compo- 
nents to the optical light curve strongly suppresses the low- 
frequency power (as they vary in anti -phase) and increases the 
high-frequency power making a bump at frequencies ^0.1-1 
Hz in the o ptical power sp ectrum (similar to that found in GX 
339^ by Gandhi etal.ll201Q) . The excess of power at higher 
frequencies makes the optical ACF narrower than that of the 
X-rays. 

The simulated CCFs for separate components and their 
joint contribution are shown in Figure ISb). The width of 
the precognition dip depends on the shape of the X-ray PSD. 
The amplitudes of both negative and positive peaks are re- 
duced, when there are two components in the optical band, 
as they enter Equation (|4]i with different signs and thus par- 
tially cancel each other. The fast rise of the CCF at zero lag 



is related to the interplay of the two optical components. The 
shape of the CCF resembles th ose found in Swift J1753.5- 
0127 ( Durant et al.ll2()()ll20TTl) . 

The corresponding optical/X-ray phase and time lags are 
shown in Figure|2fc) and (d). For the assumed disk response 
function, the phase lag c an be described by the analytical ex- 
pression (Equation 29 in lPoutan"enll2002l) : 



tanA0(/) 



rds(sina;i + X2 cosxi) 
— 1 — X2 + rds(cosxi — X2 sinxi) 



(10) 



If synchrotron dominates the optical emission (i.e., r^g < 1), 
then ^4>{f) ~ TT and At{f) oc /^^ (dashed l ines). Such a 
depen dence was observed in XTE Jl 118-1-480 jMalzac et al.l 
l2003h . In a general case, the synchrotron dominates at high 
frequencies, while at low frequencies we get 



A</>(/) 



-{Xl + X2), if Tds > 1, 



""'^d^^ f ^ ^ -f .^ ^^^^ 
■K - {XY+X2), ifrds<l. 



1 - ^-ds 



In the first case (rds > 1), this translates into constant time 
lags At « (ri + T2)rds/(?'ds - 1)- 

If the X-ray PSD is fixed, three parameters control the shape 
of the CCF and the optical PSD: r\, T2, and rds- The first two 
parameters do not affect the CCF dramatically, therefore we 
further study the role of rds- An increase of rds obviously re- 
sults in suppression of the precognition dip in the CCF (Fig- 
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ure lUb)), in the disk-dominating regime r^s ^ 1 giving the 
CCF of simple reprocessing. At the same time, the role of 
this parameter in the PSD is not so straightforward, as can be 
seen in FigureOa). The low-frequency tail is maximally sup- 
pressed at r^s ~ 1.0, at which the disk and synchrotron nearly 
cancel each other The high-frequency part is dominated by 
the synchrotron PSD and is almost independent of rds- We 
note that for rds = 1 ± a (with \a\ < 1) the PSD shapes are 
nearly the same, while the corresponding CCFs are different 
(compare the dashed and the solid lines in Figure |3ja) and 
(b)). The phase- or time-lag spectra (FigureQc) and (d)) also 
strongly depend on rds- The larger the value of rds, the closer 
the phase spectrum is to that of simple reprocessing. 



4. SUMMARY AND DISCUSSION 

In this Letter, we present a spectral model capable of ex- 
plaining optical/X-ray timing features observed in LMXBs. 
We argue that the optical emission is partially produced in the 
hot accretion flow (at distance R > 30i?s) by synchrotron 
radiation, which also provides seed photons for Comptoniza- 
tion. Additional contribution to the optical band comes from 
reprocessing of the X-ray emission in the cold accretion disk. 
We show that the presence of the two components can explain 
the observed shape of the CCF with the pronounced optical 
precognition dip. The model also reproduces another interest- 
ing feature observed in the LMXBs, namely, the optical ACF 
being narrower than the X-ray one. The explanation of the 
feature comes from the optical PSD, where the low-frequency 
part is suppressed and a bump appears at frequencies ~0. 1 — 1 
Hz, both due to the presence of the cross term. The power 
at high frequencies originates from the interplay of the two 
optical components, which also results in the fast rise of the 
CCF at zero lag. The Fourier frequency-dependent time lags 
resemble the observed ones. 

We note that the proposed model is simplified, as it does 
not account for the fact that the optical synchrotron radia- 
tion comes from the outer parts of the hot flow, while the 
X-rays are likely to originate from the very vicinity of the 
compact object. This consideration leads to two complica- 
tions: first, the synchrotron should lead the X-rays by the 
viscous timescale of about At ~ 0.1 s; second, it should 
have less power at high frequencies comparing to the X-rays. 
The latter is rather important as it leads to suppression of the 
high-frequency optical PSD, thus broadening the precogni- 
tion dip. The correct shape of the X-ray PSD is essential for 
comparison to the observed CCF, as it influences the shape 
of the precognition dip dramatically, at the same time affect- 
ing the width of the optical peak at positive lags. Generally, 
the broader is the X-ray ACF, the b roader are the two feature s 
in the CCF (as indeed is observed. iDurant et al.ll2008l 1201 ll) . 



The detailed comparison with the data will be a subject of our 
further investigations. 
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